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We introduce a new concept for the manipulation of fluid flow around three-dimensional bod- 
ies. Inspired by transformation optics, the concept is based on a mathematical idea of coordinate 
transformations and physically implemented with anisotropic porous media permeable to the flow 
of fluids. In two situations — for an impermeable object placed either in a free- flowing fluid or 
in a fluid-filled porous medium — we show that the object can be coated with an inhomogeneous, 
anisotropic permeable medium, such as to preserve the flow that would have existed in the absence 
of the object. The proposed fluid flow cloak eliminates downstream wake and compensates viscous 
drag, hinting us at the possibility of novel propulsion techniques. 
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Coordinate transformations, specifically, conformal 
maps, are widely used in fluid dynamics as a method of 
solving the biharmonic equation describing incompress- 
ible Stokes flow of viscous fluids [1 j, in domains with com- 
plicated boundary shapes. They are also applied [2] to 
solving the Laplace equation describing fluid pressure dis- 
tributions in porous media (Darcy's law), and the Poisson 
equation of electrostatics. As a subset of the separation 
of variables method, conformal transformations thus offer 
a unique and powerful approach to the forward problems 
of incompressible flow, as well as electro- and magneto- 
statics. 

The utility of more general coordinate transformations 
in the inverse electromagnetic [3H5] and acoustic [6H2] 
problems has been demonstrated recently, most impres- 
sively by showing the possibility of electromagnetic in- 
visibility, dubbed cloaking [3j HJ flQHT2] . The key idea 
that enabled the progress in those areas was the com- 
bination of coordinate transformations with coordinate- 
dependent, and often extremely exotic [5l [T3UT5] . mate- 
rial properties. Transformation optics [3H3 IIU US] and 
transformation acoustics [6-8 offer solutions to some in- 
verse scattering problems [17 by reducing them to an 
elementary one - for example, scattering off a point ob- 
ject in free space. 

There is no reason why this conceptual approach can- 
not be used in other areas of physics [18], and, in fact, it 
has already been applied to conductive heat transfer [19] , 
linear elastodynamics [5j [20H24] , surface wave [25] and 
quantum- mechanical matter wave [26] dynamics. The 
fundamental requirement for the applicability of this con- 
cept is the presence of a medium with sufficiently flexible 
properties, which enables manipulation of the coefficients 
in the equations describing the dynamical process. Here, 
we propose porous media as a substrate for transforma- 
tion fluid dynamics, and investigate the feasibility of con- 
trolling certain features of incompressible fluid flow. 

The electromagnetic cloak of invisibility, which in- 
spired this approach, manipulates the field lines of elec- 
tric and magnetic fields, and the streamlines of momen- 
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FIG. 1. (color online), (a) Schematic of the proposed fluid 
flow cloaking structure, and the definition of Problems A 
and B. (b) Solution to Problem A: pressure distribution and 
streamlines of flow, calculated using the Brinkman-Stokes 
equation ^ with Ko/b 2 — 10 -3 in Region I. The solution is 
virtually indistinguishable from the one obtained with Darcy's 
law (not shown). 



turn flux, in such a manner that these lines are virtually 
indistinguishable from those in the absence of any ob- 
ject [3]. Because the streamlines are unperturbed in the 
free space outside the cloak, electromagnetic radiation 
avoids scattering off the structure and therefore exerts 
no electromagnetic pressure on it. Achieving a similar 
level of control over the streamlines of the fluid flow out- 
side the structure would imply canceling the viscous drag 
exerted on the structure. In this paper, we demonstrate 
the possibility of fluid flow cloaking using porous media 
whose properties are more general than those considered 
in earlier studies 1271128]. 

To allow a simple analytical treatment, we consider 
a solid spherical object of radius a impermeable to the 
fluid, surrounded by a concentric permeable porous shell 
of exterior radius 6, immersed in a fluid occupying an 
unbounded region (which we refer to as free space), and 
subject to a stationary, uniform flow with asymptotic 
velocity uq in the z-direction (Fig. pTa)). For further 
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simplification, consider the flow with a small Reynolds 
number, Re = p^u^b/po <C 1, where po and po are the 
density and the dynamic viscosity of the fluid. In this 
regime, the nonlinear inertial term (u-V)u in the Navier- 
Stokes equations can be neglected, and the stationary 
incompressible flow is described by the Stokes equation: 

p V 2 u = Vp. (1) 

The fluid velocity field u is divergence- free, 

V ■ ?i = 0. (2) 

In the porous domain, the flow is described by the 
Brinkman-Stokes equation [2]: 

pV 2 u = Vp + /io^ -1 ^ (3) 

where p is the effective viscosity of the flow in the porous 
medium, and n is called the permeability. 

In general, for an anisotropic porous medium, ft -1 is a 
rank- two tensor [2], and p is a tensor of rank four. While 
the implementation of anisotropic permeability using lay- 
ered media is well-known (see §5.8.3, "A Layered Medium 
as an Equivalent Anisotropic Medium" in Ref. [2]), effec- 
tive viscosity of fluids in porous media has been a subject 
of some debate [28]. To evaluate the degree of control 
that can be achieved by manipulating only permeabil- 
ity, in this study we assume that the effective viscos- 
ity in the porous medium (p) is isotropic, constant, and 
equal to the fluid viscosity in the absence of a porous 
medium (/io); this assumption is consistent with estab- 
lished practice [28]. Whenever a numerical calculation 
occurs, we assume for concreteness that po = 10 3 kg/m 3 
and /io = 10 -3 Pas, which approximates water at room 
temperature. Additionally, it is assumed that the radius 
of the structure b = 1 mm. 

In the limit of small permeability, n <C 6 2 , the term 
/iV u in Eq. (|J) can be neglected in comparison with the 
term p^n^u, and the Brinkman-Stokes equation reduces 
to Darcy's law [2], 

Vp -poK^u, (4) 

which, combined with the continuity equation for an 
incompressible fluid gives Darcy's pressure equa- 
tion 0: 

- VftVp = 0. (5) 

If we allow the porous medium to fill the entire space 
around the solid sphere, Darcy's equation ([5| is at least 
approximately valid in all regions of the flow. We then 
formulate Problem A as follows: assuming that the 
porous medium in Region I (r > b) is uniform and 
isotropic with permeability = const, and that the per- 
meability tensor n is uniaxially- anisotropic and radially- 
symmetric in Region II (a < r < 6), i.e. 

n = n r (r)rr + n t (r)00 + Kt{r)(jxj), (6) 



find the functions n r (r) and K t (r) such that a solution 
exists in which the velocity is constant everywhere in 
Region I: 

u(r, 0) = uoz = uo(r cos — sin#). (7) 

Note that because Region I contains a porous medium, 
there is a pressure gradient in the ^-direction even for 
a uniform ("plug") flow ([7|, in which the viscous fluid 
stress is absent. 

The formulation of the more practically important 
Problem B is the same as that of Problem A, except 
that Region I (r > b) contains no porous medium. In 
both Problems, we require the standard, no-slip bound- 
ary condition, that is, u = at r = a, which corresponds 
to a surface of an impermeable solid (Region III). 

Because Problem A in Darcy's approximation is math- 
ematically equivalent to the problem of cloaking a spher- 
ical object from a uniform electrostatic field (except for 
the boundary condition), we could simply use the solu- 
tion obtained in Ref. [3], 

/ \ b f r — a\ 2 . x b 

Kr{r) = ^0 7 , ^W=/^07 , (8) 

b — a \ r J b — a 

which corresponds to the following transformation: 

r' = b{r - a)/(b - a), 6' = 0, <jf = 0. (9) 

While the boundary condition on the surface of the 
cloaked object (r = a) used in electrostatics is different 
from the no-slip boundary condition for the fluid flow, in 
three dimensions a spherical cloak described by Eq. ([8| 
remains a perfect cloak regardless of the kind of homoge- 
neous boundary condition used at r = a; see Ref. [12] for 
a discussion of the boundary conditions in two and three 
dimensions. 

The pressure distribution and the streamlines of plug 
flow around the structure obtained in the approxi- 
mation of Darcy's law are shown in Fig. [TJb) , proving 
that the distributions ([8| represent an exact solution to 
Problem A. The plotted solution was calculated through 
the axisymmetric formulation of the Brinkman-Stokes 
equation available as the porous flow model in the finite 
element analysis package COMSOL Multiphysics [29] . 
Anisotropic permeability in the Brinkman model was im- 
plemented by customizing the stress tensor prescribed by 
the standard model. When the permeability is small 
enough, the viscous contribution to the fluid stress is neg- 
ligible compared to the Darcy's force, and the solution of 
Brinkman-Stokes equation is indistinguishable from the 
solution of the Darcy's equation. 

Now, we turn to the more interesting Problem B where 
the unbounded Region I is free of porous medium (see 
Fig. []Ja)), and the flow is described by the Stokes equa- 
tion in Region I and the Brinkman-Stokes equation in 
Region II. The particular case of a creeping flow around 
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FIG. 2. (color online). Solutions to Problem B: maps of the 
axial velocity normalized to velocity at infinity (u z /uq) and 
streamlines of flow around and through the porous spherical 
shell. Left to right: Re < 0.5 (a), Re = 2.5 (b) and Re = 4.5 
(c). 



a uniform, isotropic porous shell has been studied by var- 
ious authors using the Stokes stream function [27JI2H] and 
other methods [30] [31]. Here, we generalize the known so- 
lutions [271 HE] of the forward problem to allow the inclu- 
sion of inhomogeneous, although still radially-symmetric, 
uniaxially anisotropic permeability (|6|, and use them to 
solve the inverse problem (Problem B). 

In Region I (absent of any porous medium), the most 
general solution to Eq. (|3| that satisfies the condition of 
plug flow at r — > oo can be written as [28 



u r = Uq cos ( 
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where A and B are some coefficients that depend on the 
details of the flow in the porous Region II. In the exact 
cloaking solution we are seeking, both A = B = 0; this 
means that all components of the viscous stress tensor, 



djUi) 



as well as the hydrostatic pres- 



sure p (gauged to p = at infinity) vanish everywhere in 



Region I, leading to zero fluid stress 



-pfii: 



and a vanishing drag force. The latter can be calculated 
as F z = § Sr (j Z jdSj, where Sr is a sphere of radius R. It 
is easy to see that in the limit R — >• oo only the terms 



~ B/r in Eq. (10) give finite contributions. Thus, the 
drag force is proportional to the coefficient B and is en- 
tirely independent of A [28]. Therefore, for the purpose 
of eliminating only the drag force, but not the laminar 
wake (deviation of the flow from uniformity), a weaker 
form of the cloaking problem can be postulated that only 
requires B = without a constraint on A. This wider 
class of "weak cloaking" solutions is not analyzed in this 
work. 

To find the exact cloaking solution (A = B = 0), we 
employ numerical optimization of the Brinkman- Stokes 
equation ([3| with the unknown anisotropic permeabil- 



ity profile using COMSOL [29]. The inertia terms in 
both the Brinkman and free-space Navier-Stokes equa- 
tions were neglected. The optimization solver discretizes 
the unknown one-dimensional functions of the spherical 
radius, ^~ 1 (r) and ^^ 1 (r), and treats the nodal values of 
these functions as the optimization variables to solve for. 
The optimization goal minimized by the solver is chosen 
as A 2 + B 2 \ the fully converged solution has A « B « 
with a numerical accuracy of better than 10 -8 . Conse- 
quently, the flow pattern everywhere outside the porous 
shell satisfies the requirement u = const, as seen from 
Fig. Ufa). 

The distribution of k components necessary to imple- 
ment this flow pattern using a shell with the aspect ratio 
b/a = 2 is shown as the thick solid line in Fig. [3] The 
function n r (r) is negative in its entire range a < r < 6, 
and the function n t (r) is only positive for ro < r < b and 
negative for a < r < ro, where ro ~ 1.8a. The medium 
with a negative permeability is an active medium: it must 
use external energy to provide acceleration rather than 
deceleration to the fluid permeating through the porous 
medium. A physical implementation of such a medium 
can be foreseen as a distribution of minute pumps that 
help propel the fluid and compensate in a controlled fash- 
ion the viscous pressure drop that exists due to the ve- 
locity gradient in the a < r < b region. Potential can- 
didate technologies for such a directional acceleration on 
the micro-scale level are electro-osmotic [32], pneumatic 
and piezoelectric [33] micro-pumps. 

So far we have been neglecting the nonlinear advective 
term in the Brinkman equation, which is only valid for 
Reynolds numbers i?e « 1. In this regime, the perme- 
ability profile is independent of the asymptotic velocity 
Uq. Using numerical optimization of the full Brinkman 
and Navier-Stokes equations, 



Po(u • V)u = — Vp + /ioV U — juLQK 



(11) 



where n 1 = in Region I, we have also investigated 
the possibility of higher- Re flow cloaking. In Eq. (11), 



we assume that porosity in Region II, e p , equals unity, 
which is consistent with our assumption \i = /i - 

Figures [2|a-c) demonstrate that the flow can be 
cloaked for Reynolds numbers up to approximately 4.5; 
at larger Re the assumption that the flow in Region I 



is given by Eq. (10) is no longer accurate, and the op- 



timization goal based on the A and B coefficients loses 
its meaning. The drag force F z , normalized to the ideal 
Stokes drag for a sphere of radius a (F$t = ^TTfiuoa pQ), 
was estimated by direct integration of the fluid stress, 
and shown to be < 0.01 at all Re < 1, w 0.157 at 
Re = 2.5 (Fig.^b)), and w 0.48 at Re = 4.5 (Fig. [^c)). 
The resulting permeability profiles for various Re values 
are shown in Figj3j We observe that these profiles are 
virtually independent of the velocity magnitude up to 
Re ~ 2.5, but above that value they evolve quickly with 
Re. 
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FIG. 3. (color online). Solutions to Problem B: radial (top) 
and azimuthal (bottom) components of inverse permeability 
(normalized by = a 2 /4) as functions of the spherical radius 
r corresponding to solutions of the exact cloaking problem 
(A = B = 0) with various Reynolds numbers. 



In conclusion, we have postulated and solved the prob- 
lem of cloaking a spherical object from the fluid flow in a 
uniform porous medium (A), and in free space (B), using 
a permeable porous shell. The new concept of utilizing 
permeable media for fluid flow management opens the 
door to novel hydrodynamical approaches to compensat- 
ing the viscous drag force, eliminating laminar wake and 
inhibiting the formation of turbulent wake past moving 
bodies. The emergence of these strategies hints at the 
possibility of wake- free distributed propulsion systems. 
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